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Abstract. Activation and deactivation by inelastic collisions have been extensively studied at 
unimolecular reactions in gas phase where they are crucial for equilibration. As equilibration 
means an increase of entropy, the mechanism can also be considered responsible for entropy 
production. The numerous results of theoretical treatments show a remarkable agreement 
with experimental data. Analysis of the statistical physicochemical argumentation reveals the 
stochasticity of quantum transitions as the essential assumption. Under this premise, master 
equations have been used that are known to deliver equilibria and entropy production. 

Here we examine the hypothesis that the ubiquitous inelastic interactions in gas and liquid 
phase may represent a source of entropy beyond chemical reactions, rotational activation/deac- 
tivation being the prevailing mechanism in gas dynamics at room temperature. For a quantum 
mechanical two-state model the master equations are formulated which yield entropy produc- 
tion and equilibration in translational degrees of freedom until, at conserved energy, a station- 
ary Maxwell-Boltzmann distribution is reached. 

The relaxation rates show features that can be checked by monitoring thermal relaxation 
in gas phase. Depending on the composition, first or second order processes are predicted. 
The temperature dependence of relaxation rates is determined by the activation energy of the 
dominating quantum transition. Thus, experimental verification will allow to decide to which 
extent this hypothesis describes thermal relaxation. It would support a connection between the 
macroscopic second law of thermodynamics and the microscopic stochastic collapse in quan- 
tum mechanics. They are both experimentally secured facts which in theory emerge as special 
elements beyond Hamiltonian dynamics. 

It is also shown that inelastic collisions are connected with velocity-dependent forces and 
result in a new analog of the Fokker-Planck equation where they replace Langevin dynamics 
as a non-Hamiltonian dissipative mechanism. 

Keywords: Entropy production, dissipation, Fokker-Planck equation, master equation, inelas- 
tic collisions, gas dynamics 
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1. Introduction 

Recently, several new attempts were made to prove second-law conform behavior of classical 
systems when work is applied [HI El El or at sudden release of a constraint [4J. It is in 
fact possible to derive an increase of entropy in such situations using classical mechanical 
statistics. On the other hand, it is known that Gibbs entropy - which is always referred 
to here - is time-invariant under classical Hamiltonian dynamics and cannot change in the 
situations mentioned. The seeming contradiction is resolved by the observation that the 
above theories are all based on the assumption of equilibria, which are known to be states 
of maximum entropy Q. For instance, in our own approach it is proven that immediately 
after release of a constraint, any system is found in a state of higher entropy - provided it is in 
equilibrium at the start and the new position [4J. However, when it is assumed that a process 
ends with equilibrium, the treatment tacitly departs from the field of Hamiltonian dynamics 
without declaring the new theoretical basis. One has to conclude that the cited approaches are 
formally correct, but the physical mechanism of entropy production during equilibration is 
still the ultimate unsolved problem. Moreover, it has to be described using an adequate non- 
Hamiltonian ansatz. An example of this kind is the derivation of Jarzynski's work theorem 
from Langevin dynamics, i.e. a non-Hamiltonian phenomenological ansatz J6]|. 

Because of the maximum entropy property of equilibria, relaxation mechanisms are 
candidates for the explanation of entropy production. Equilibration by activating and 
deactivating inelastic collisions were extensively studied at unimolecular reactions in gas 
phase B71[8l|9l. These are reactions of the type 

M + A^±M + A*^M + P (1) 

where molecules A are activated/deactivated by collisions with molecules M, and can decay 
from the activated state A* to product(s) P. A and M may be the same molecular species. The 
last step is responsible for the denomination as a unimolecular reaction. It requires 'inter- 
nal energy randomization' rTTOll for transferring energy to the crucial reaction coordinate. The 
results of numerous theoretical treatments show a remarkable agreement with experimental 
data. Analysis of the statistical physicochemical argumentation reveals the stochasticity of 
quantum transitions as the essential assumption. Under this premise, master equations have 
been used that are known to deliver equilibria and entropy production. Rotational and vi- 
brational activation were considered for the step A ^ A* with experimentally determined 
inelastic cross sections [9|. The final decay also needs deactivating collisions as indicated 
by pressure dependence of the transmission coefficient of the decay which can increase with 
pressure due to increasing frequency of collisions [fTT||. 

The successful description of inelastic collisions as stochastic events by master equations 
can be regarded as an approximation that requires justification by invoking decoherence phe- 
nomena lfT2lllfT3l . Here we take stochasticity for granted. We reconsider the old experiments 
and theories as they suggest a very general role of such collisions. The proper hypothesis 
of this work claims that the ubiquitous inelastic collisions in gases and liquids provide an 
entropy-producing mechanism already under the threshold of chemical reactions. As will be 



2 



Figure 1. The principle of inelastic QM/MM interaction. At collision of molecules, energy 
transfer between translational and internal degrees of freedom by inelastic processes can occur, 
which are the focus of this investigation. 

shown, there is a realistic way for experimental verification. 

The hypothesis of entropy-producing inelastic collisions is substantiated by the existence 
of QM energy levels that can be excited even at low and medium temperatures. The proper 
verification/falsification of the hypothesis requires prediction of measurable features. This is 
done by theoretical treatment of an internal two-state model which is shown schematically in 
Fig. [T] It possesses the essential features of relaxation mechanisms employed elsewhere IHIIl 
and allows to derive relaxation rates. 

Theories of unimolecular reactions focused on the equilibration among molecular species, 
i.e. chemical relaxation. Our hypothesis concerns thermal relaxation of both internal and 
translational degrees of freedom of the colliding molecules, which were not considered at 
unimolecular reactions. Translations are denoted and treated as classical (or molecular me- 
chanical 'MM') while the QM degree of freedom are represented by two distinct energy levels. 
We investigate the impact of the stochastic QM transitions on the MM system of translations. 
There are also theories which take into account the influence of a stochastic bath on a QM sys- 
tem, e.g. by applying a stochastic Liouville equation IIT4"l[T5Tl to spins or considering quantum 
information [ 16 1. In this study, however, the QM system is the 'stochastic bath' for classically 
moving particles. 

The article is organized as follows: In section 2, we examine available inelastic processes 
which are neglected or inadequately treated in usual classical dynamics. The master equa- 
tions of the two-state model are developed and transformed to equations for various entropies 
in section 3 where also rate constants are derived. In the following section 4 we propose in- 
corporating the inelastic process in the Fokker- Planck equations as a dissipative mechanism 
which can replace the phenomenological Langevin dynamics. In the final discussion, the fo- 



cus is placed on the justification of the underlying assumptions and formalism, and the relation 
to Shannon information. 



2. Inelastic processes 



When two neutral molecules collide without chemically reacting as shown in FigJTJ part of 



the kinetic energy E\ + E 2 of the centers of mass is transiently transformed into internal 
energy E int , for instance, by induced polarization. Here, internal energy denotes all kinds 
but translational energy. The resulting van der Waals interaction is the only effect at slow 
collisions below the threshold of excitations. Persistent exchange of energy with rotational or 
vibrational levels takes place when molecules are in excited states or kinetic energy is large 
enough to allow excitation. A quantum mechanical calculation proves that internal energy 
remains unchanged only when the collision is very fast, the so-called rapid-passage case [17J. 

Here, internal energy comprises all quantized degrees of freedom with discrete spectra 
including rotations. At room temperature 300 K corresponding to a wave number 
v(k B T) « 200cm 1 , rotations (v R < 20cm' 1 ) play a major role in gases and liquids. 
In liquids, also librations (v L « 300 - 900cm' 1 in H 2 at 0C°) and hindered translations 
(v L « 200cm" 1 in H 2 0) contribute to energy exchange at collision. Monatomic and most 
homonuclear diatomic molecules are rotationally inactive and provide no suitable internal 
degrees of freedom at 300 K. Here, only impurities and collisions with walls confining the 
gas present mechanisms for inelastic energy exchange. Moreover, vibrations are available 
with low frequencies at larger molecules, and at higher temperature electronic transitions are 
induced everywhere by collisions. The variety of low energy transitions available in gases 
and liquids that can be excited at room temperature demonstrates the importance of inelastic 
interaction in the classical regime. As a consequence, a finite fraction of interacting particles 
will, after a while, be found in some excited state, and a permanent uptake and release of 
energy by inelastic processes is taking place. This side-effect of intermolecular interaction 
and its implications are the subject of this study. 

It is worth considering the reason why entropy is invariant when the dynamics is governed 
by a Hamiltonian H(q,p) = K(q,p) + V(q), (q,p) being positions and momenta. For the 
classical Gibbs entropy 
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kg ist the Boltzmann constant and P the probability density. The second form ([4]) obtained by 
partial integration proves that only a momentum-dependent force p = F(q, p) could influence 
entropy, but not a force derived from a potential function V(q) as assumed in Hamiltonian 
dynamics. There also the energy H(q,p) is invariant. On the other hand, inelastic interactions 
stand for uptake or release of energy portions AE in short time intervals At. Because kinetic 
energy changes according to dE = vdp = vFdt, one finds that during the transition a force F 

1 AE 

v At 

is acting. The significance of this relation lies in the fact that here the force indeed depends 
on velocity v or momentum p. This is a strong indication for entropy production by inelastic 
processes that can be extended to a mathematical argument. Here, however, the final proof 
will be given by a derivation that starts from ([3]). Velocity-dependent forces are characteristic 
of Langevin dynamics. In simulations, they occur in thermostat algorithms where they mimic 
heat exchange with an environment. 



3. Master Equations of a two-state model 

We now consider a two-state model according to Fig|T]The following assumptions are made: 
(i) Atoms and molecules are point particles with classical probability density functions (PDF) 
P (1) (x, p, t) for a single particle and P (JV) (x,p, t) for an N-particle system together with the 
corresponding Gibbs entropy 0. (x, p) are the respective Cartesian coordinates and momenta 
of the centers of mass. A classical treatment is justified by the uncertainty principle which 
predicts classical behavior for a box length Ax/nm » (M/(g/mol))~ l , i.e. already in the 
higher nanometer range at any molecular mass M. (ii) In the classical world QM transitions 
can be taken into account as forces acting for a period of time or by discrete energy portions 
when they couple to translation. As a matter of fact, molecules and single atoms in ordinary 
experiments at room temperature are at the interface between the classical MM and QM 
regime: their rotational or electronic degrees of freedom require a QM treatment, (iii) The 
transitions are random, stochastic processes, which justifies a master-equation approach, and 
(iv) energy exchange is assumed to occur strictly without net change in energy on the average 
and therefore without transfer of entropy. 

We neglect the simultaneous energy transfer on both colliding particles. Instead, we consider 
the action of a generic, non-specified QM system on single particles that feeds in energy 
packets (power strokes) +AE > or removes energy in portions -AE. The symbol AE = \AE\ 
stands for the absolute value, E = Ekm = p 2 /2m is the kinetic energy of translation. As 
the events of energy exchange are random, stochastic QM transitions we employ a master 
equation for the classical probability density function to account for their action on translation. 
Elastic interactions are not included from the beginning, but added later. We start with the 
development of the energy density P(E, t), < E < oo, which facilitates the connection with 
the QM system. 

Let us assume that the uptake of AE occurs with rate w + and the release of the same amount 
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of energy with rate w . Uptake of energy simply means that in a time interval dt, part of the 
density is lost everywhere and shifted to higher energy, where it appears at some E > AE, 

dP upt (E, t) = (-w + P(E, t) + w + P(E - AE, t)) dt if E > AE 

dP upt (E, t) = -w + P(E, t) dt if E < AE (6) 

A release of -AE takes place only if the energy is available, i.e. for E > AE. Therefore, one 
has two cases 

dPrei(E, t) = (-w~P(E, t) + w~P(E + AE, t)) dt if E > AE 

dP rel (E,t)= +w~P(E + AE,t)dt if E < AE (7) 

In any case, part of the density is shifted to lower energy during release, see also Fig. [2} 

By integration over energy it is easily seen that the normalization of probability density 
is always maintained. Taking both processes together and making use of the Heaviside step 
function = Q(E - AE), one obtains 

dP(E, t) = [-(w + + ®w )P (E, t) + &w + P (E - AE, t) 

+ w~P(E + AE,t)]dt (8> 



Because of the different behavior at low energy in (|6jj7j), the rates w + and w cannot be equal, 
but differ by a factor ip = w~/w + . It is assumed that there is no net change of energy, i.e. 



oo oo 

J P(E,t)dE = w + J P{E,t)dE 



w~ P{E,t)dE = w + P{E,t)dE (9) 



AE 



which implies w + < w and (p = J P(E,t)dE. This condition will be employed in the 

AE 

appendix for calculating the time dependency of entropy at constant energy. 

It is obvious that this dynamics is suited to broaden an initially narrow distribution or 
constrict a broad distribution, and possibly produce entropy. This is rigorously proven in the 
appendix for the time derivative of the entropy S E belonging to P(E, t), 

oo 

= ******** do) 



and, more importantly, for one-particle phase-space Gibbs entropy defined by 

oo 

S (l) = -k B J P m ( V ,t)\nP w (v,t)d 3 p (11) 

o 

of the one-particle phase-space PDF P (1) (p, t) on the assumption of energy conservation. The 
result is dS^/dt > 0, i.e. entropy never decreases. More precisely, entropy turns out to 
increase as long as the stationary state is not reached where entropy remains constant. In 
other words, coupling to and energy exchange with a quantum system is causing dissipation 
and equilibration in a classical system. 
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Figure 2. Broadening of probability density by release (left) and uptake (right) of energy AE, 
Fractions of the PDF are shifted to the flanks whereas the original PDF (bold) is flattening. 
Note that systems with E < AE cannot give away their energy. They are not damped, but can 
be speeded up. 



The stationary solution is obtained in the same way (see appendix) or by solving eq. 
(|6jj7]) with the condition dP = dP upt + dP re i = 0. The unique integrable, real solution of 
either approach is a decreasing exponential which is actually independent of AE although, 
according to Q7]), the two intervals have to be treated separately. It is also independent of the 
transition rates w +/ ~, which confirms its equilibrium character. A posteriori, one can write it 
as a Maxwell-Boltzmann distribution 

P[ 1) (E)=/3cxp(-/3E) (12) 

with /3 = 1 /k B T which is characterized by mean energy or temperature T. This finally proves 
that stochastic interaction with a QM system is an entropy-producing physical mechanism 
that generates the well-known equilibrium distribution Ps(E) of maximum entropy. This 
also holds for the one-particle PDF P^ (/?), see appendix. The mechanism proposed here is 
thus equivalent to the phenomenological Langevin mechanism that is usually assumed to be 
present and actually causes convergence towards equilibrium. 



According to eq. A. 6 entropy production and equilibration occur with a rate w + = ipw , w 



being the rate of energy transfer into the quantum system. Energy conservation ([9]) together 



with ( 12 ) reveals <p to be a Boltzmann factor <p = exp (-fiAE), which indicates also equilibrium 
in the QM system. This finally allows application to inelastic collisions A + A ^± A + A*. 
Collision theory predicts a pressure dependence of second order and hence a relaxation rate 

w + oco- in ex V (-J3AE)p 2 A (13) 

for a homogeneous gas with pressure p A if the molecules can be excited/deexcited. If there is 
a major rotationally inactive component / with pressure pi that offers no quantum transitions 
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at a given temperature, one expects 

w + oc a in exp (-(3AE) p A pi (1 4) 

where p A is the pressure of an active impurity and cr in the respective inelastic cross section. 
These relaxation rates possess characteristic properties that can be checked experimentally, 
namely a typical temperature dependence given by the threshold energy AE of the (lowest) 
quantum transition and a pressure dependence due to the composition of the gas. Relaxation 
is expected to be strongly hampered at rotationally inactive gases. The two-state model will 
be valid at sufficiently low temperature where higher states are not excited; it will be ap- 
proximately valid where such processes take place and a larger relaxation matrix is required. 
Resonance effects will not alter the essential features of the relaxation rate. It was tacitly as- 
sumed that the above inelastic cross section is an average, and likewise that a temperature can 
be assigned to the system during relaxation. 



4. Analog of the Fokker- Planck equation 

The full time-dependence of the PDF P (1) (x, p, t) is usually described by a Fokker-Planck 
equation (FPE), assuming that dissipation is given by the Langevin equation of motion, 
x = -^x + A(t), with friction constant $ and a stochastic force A(t) with vanishing mean. 
The generalized FPE or Chandrasekhar equation [18] is a differential equation of the form 

p(i) = Rp(i)+DP (1) (15) 

containing two differential operators, R = v • V v + (F(x)/m) • V v and D = gV v (\ + (k B T/m)V v 
with v = x. RP (1) is the regular part following from Hamiltonian dynamics with force 
F(x) = -V x V(x) derived from an external potential. The term DP being proportional to the 
friction constant g, is non-Hamiltonian and accounts for dissipation according to the Langevin 
equation. 

In contrast, we here investigate dissipation basing on inelastic interaction with a QM 
system. Equation ^ can be given the shape of an integral operation writing P = I E P. Then 
the respective change of P (1) can be written as P { q M / MM = I P (l \ where I is the integral operator 

oo 

IP (1) (x, p, t) = - f [-(w + + ®w~)8{E' -E) + @w + 6(E' - E + AE) 

PJ (16) 

o 

+ w'6(E' -E- AE)]E'P (i \x, p', t) dp' 

Here E = p 2 /2m and E' = (p') 2 /2m. The direction of momentum is not affected. The 
transition from P to P (1) is derived in the appendix. The operator I E , and hence also I, do 
produce entropy and have a stationary solution which coincides with the Maxwell-Boltzmann 
distribution. It depends on the rates w +/ ~ and threshold energy AE and can replace the 
Langevin operator D to account for dissipation. Therefore, we propose a new equation that 
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is formally analog to the generalized FPE P (l) = RP' 1 ' + D P w and describes the dissipative 
part of particle interaction, namely 

= Rp(D+ip(i) (17) 
It has the same regular part KP^\ but a different dissipative part that is typical of 



inelastic QM/MM interaction. Note that ( 17 ) does not fully replace the FPE equation because 
elastic interactions are excluded for the time being. The stationary solution of the new 
equation is 

4 1} (x,p) = Ns l exp [-/?(p 2 /2m + y (x ))] (18) 

like for the FPE. TVs is a normalization constant. The stationarity with I was proven in the 
previous section, and is easily seen for R by using the above explicit form of the operator. 

The new integro-differential equation ( [17] ) for a single particle is straightforwardly 
generalized to N-particle systems with x = (xi...x#),p = (Pi-.p^) and a density P (Af) (x,p). 
Defining the operators R, and I, as one -particle operators like above which act each only on the 
phase-space coordinates x, and p ; , one constructs the N-particle operator as a sum 2,-(R,- + I,) 
that governs the motion of the N-particle density as 

P m = ^(R,- + IdP (N) (19) 

i 

The increase of entropy due to the integral terms I, follows from the behavior proven for the 
one-particle density , while the regular parts R, do not contribute. The stationary solution is 
the Maxwell-Boltzmann distribution 

3N 

P ( s N) (x,v) = P s (x)Y\py(V i ) 

i=l 

P ( S(Pi) =N- 1 eKv(-J3pt/2m i ) 

/f°(x) = jV 1 exp(-/?V(x)) (20) 

which is known to be the maximum entropy PDF. The /V, are normalization factors. Here, V(x) 
also includes all elastic interactions between particles and external potentials. They were not 
yet included in the one-particle equation ( [T7j ), but occur only in the N-particle system ( 19 ). 
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5. Summary and discussion 



We have formulated the hypothesis that the ubiquitous inelastic collisions in gas and liquid 
phase represent a source of entropy. Stochasticity of the underlying quantum transitions 
is the essential assumption to be made. It is supported by the numerous applications to 
unimolecular reactions where the experiments are well described. Rotational levels are 
available for excitation at ordinary temperature, followed by vibrations. We have proposed 
to check the hypothesis experimentally by measuring (or evaluating existing data about) 
thermal relaxation in gas dynamics for which the variation with temperature, pressure and 
composition ( 13|14 ) was derived from a model. Note that the dependence on pressure and 
composition, for instance the presence of inert gases, was already analyzed at unimolecular 
reactions (5]|. Apparently, it is possible to answer the important question of entropy production 
experimentally. The optimum temperature range should be low enough (a) to avoid chemical 
reactions and (b) to facilitate evaluation by restricting excitation to essentially one single 
energy level. 

It is often believed that elastic collisions of molecules in a gas provide a mechanism for 
thermal relaxation, i.e. emergence of a Maxwell-Boltzmann distribution independent of the 
starting conditions. The Boltzmann equation is occasionally invoked as an approximation 
for dilute gases, and generalized for dense systems lfT9ll . In principle, dephasing cannot be 
excluded to yield an entropy-producing master equation as part of the equation of motion of 
the density matrix as shown for model systems lfT3l . On the other hand, quantum mechanical 
potential scattering is more closely related to its classical counterpart than inelastic scattering, 
and classical potential scattering is entropy-conserving in arbitrary systems. In the end, the 
role of elastic scattering is also a question that can be answered by the proposed experiments. 

If there is a stochastic process which interferes with Hamiltonian dynamics and can 
be described using master equations, then it is known to imply entropy production for 
mathematical reasons. This can be understood by invoking the effective identity of entropy 
and (missing or needed) Shannon information. Stochastic events hamper the prediction of the 
future of a system from its present starting conditions. Each time, extra information has to be 
added for a full description of the system developing in time. In a discrete representation, 
the string that fully encodes the system in phase space thus becomes longer and longer. 
This holds similarly for the PDF with the reservation that information has an upper limit 
dependent on energy. The second law of thermodynamics clearly states that actually this 
happens permanently in reality. 

A common way of introducing dissipation is the use of phenomenological Langevin 
dynamics which adds stochastic and velocity-dependent forces to Hamiltonian dynamics. We 
have shown in passing that in fact stochastic energy transfer by inelastic collisions results 
in velocity-dependent forces. The corresponding terms can fully replace the Langevin terms 
in the Fokker-Planck equations to induce dissipation and relaxation to Maxwell-Boltzmann 
distributions. This has far-reaching consequences for the explanation of numerous relaxation 
phenomena. 
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Up to this point, is has become clear that - according to our hypothesis - inelastic collisions 
are candidates for causing entropy production. Moreover, this suggests that quite generally 
energy exchange of any QM system with a MM system via discrete transitions as illustrated 



question when and why the QM transitions can be treated this way, i.e. as stochastic power 
strokes for the classical system. This touches the explanation of the transition from QM to 
classical dynamics, which is out of the scope of this work. It is still under debate and subject 
to experiments iPlP . 

The available experiments on unimolecular reactions already suggest two interesting 
conclusions that could be corroborated by the proposed new experiments. Firstly, they 
indicate that colliding molecules in large systems - where the reactions are satisfactorily 
described - behave like under measurement conditions. Transitions between energy levels 
involve the collapse of the wave function occurring with a probability that is well described 
by QM, but the unpredictability of the time point of the collapse makes them a perfect random 
generator. Secondly, there seems to exist a direct connection between macroscopic entropy 
production as stated by the second law of thermodynamics and the microscopic stochastic 
collapse in quantum mechanics. Apparently, stochasticity can explain the increase of Shannon 
information or physical entropy. This would connect two experimentally secured facts, 
quantum transitions and the second law, which in theory emerge as special elements beyond 
Hamiltonian dynamics. 
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Appendix A. Time derivative of entropy 

We first derive the time derivative of the entropy S E , 



by Fig{T]will produce entropy in both systems like in the model. What remains open is the 




(A.1) 



o 



using 



P(E, t) = (-w + - ®(E - AE)w~) P(E, t) +w + Q(E - AE) P(E -AE, t) 



p p- 



(A.2) 



+ w~ P(E + AE, t) 
■* ' 
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from The one-particle entropy will be treated in a second step. As the integral over P 
was already shown to be conserved, it suffices to consider 



S E /(k B w ) 



oo 

-f 

w J 



PlnPdE 



= (pj PlnPdE + J PlnPdE -<p J P'hiPdE-J P + InPdP 

AE AE 



(A.3) 



This is where energy conservation ([9]) comes into play. In contrast to P and P , P + is not 
normalized in (0, oo). The corresponding normalized PDF is P = P + /(p because energy 

OO OO 

conservation demands J P + dE = J PdE = if. Then one finds 
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AE 



C = J P lnPdE = J PlnP + dE = J Pln(P<p)dE 

AE 

oo oo 

p 



= J PlnPdE + J PlnPdE + ln<p 



o 







Pln-dE +A + lnw 
P Y 



(A.4) 



and 



(fi(A-C) = -(p\n<p + ip J PlnPdE 

o 



CO CO CO 

J PlnPdE = J P + \nP + dE = (p J 



B = I PlnPdE 

AE 



dE = (p P\n(P<p)dE 



J PlnPdE + \n(p 

L0 



(A.5) 



D = J P + \nPdE = <f J PlnPdE 

o o 



B- 



D = <fiJ PlnPdE + ipln<p 

o 
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By adding (A. 4) and (A. 5) and using Gibbs' inequality for arbitrary PDF's, J dEPXn ~ > 

o 

one finally obtains 



S E /(k B w-) = (p(A - C) + (B - D) = ip 



CO CO 

J P\njdE + J Pln^dE 



> 



(A.6) 



This proves S E > 0, i.e. this entropy never decreases. Gibbs' inequality further implies that 
entropy is constant if and only if P = P, which is the condition for the stationary solution. 
Because of the above definition of P, this is equivalent to the condition P(E + AE) oc P(E) that 
is satisfied only by an exponential or constant function. In order to be integrable, it must be a 
decreasing exponential, which proves that the stationary solution is the Maxwell-Boltzmann 
distribution. In summary, one can state that entropy S E always increases, S E > 0, until the 
stationary state is reached. S E is proportional to w + = w~(f, which is the rate for the uptake of 
kinetic energy. 

The interesting one-particle phase-space entropy 5 (1) belongs not to energy, but the PDF 
of momentum ,P (1) (p, t) = (4nmp)~ l P(E, t). This follows from E = p 1 12m and dp = dE m/p. 
Integration over momenta confirms the correct normalization 



o o 
The connection between the two entropies is achieved by 

oo 

S (1) = -k B j P (1) (p, t) In P (1) (p, t) d 3 p 
o 

OO 

= -k B J (AnmpY l P(E, t) [lnP(£, t) - \n(Anm) - In p] Anp 2 [dE m/p] 



(A.7) 



o 



(A.8) 



= -k B | P(E,t) 
o 

1 



In P(E, t) - ln(47rm) - \ In 2mE 



dE 



= S E + -k B (In E) + const 

It seems reasonable to accept (lnis) = const as equivalent to energy conservation. With 
this assumption, one finally obtains the desired time-dependence of one -particle phase-space 
entropy S 1 , namely > 0. Also the entropy always increases, > 0, until the 
stationary Maxwell-Boltzmann distribution is reached. 
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